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A series of calculations on MgE>2 and related isoelectronic systems indicates that the layer 
of Mg 2+ ions lowers the non-bonding B n (p z ) bands relative to the bonding a (sp x p y ) bands 
compared to graphite, causing a — > n charge transfer and a band doping of 0.13 holes/cell. 
Due to their two dimensionality the a bands contribute strongly to the Fermi level density of 
states. Calculated deformation potentials of F point phonons identify the B bond stretching 
modes as dominating the electron-phonon coupling. Superconductivity driven by a band 
holes is consistent with the report of destruction of superconductivity by doping with Al. 



The dependence of the superconducting critical 
temperature T c on structure is unclear after decades 
of study, and recent discoveries have further con- 
fused the issues. In conventional superconductors 
(the high T c cuprates comprise a special case) it has 
been generally thought that high symmetry, prefer- 
ably cubic, is favorable for higher T c . This trend has 
held up in elemental superconductors (Nb, La under 
pressure, at 9 K and 13 K, respectively), for binaries 
(Nb3(Al,Ge), 23 K (!]]), (pseudo)ternaries such as 
(Ba,K)Bi03, J| and even fullerene superconductors 
(T c to 40 K §). 

Recently several intriguing counterexamples to 
this trend have come to light. Semiconducting 
HfNCl (and ZrNCl) is a van der Walls bonded set 
of covalent/ionic bonded layers, but superconducts 
up to 25 K when doped (intercalated) with alkali 
metals. j| There is evidence that the surface layer 
of solid Ceo becomes superconducting up to 52 K 
when it is injected to high hole concentrations. || 
And most recently, it is reported that the layered 
metal/metalloid compound MgB2 superconducts at 
^40 K, HI?]] which is by far the highest T c for a bi- 
nary system. The B isotope shift of T c reported by 
Bud'ko et al. and most other early experimental 
data H suggests conventional BCS strong-coupling 
s-wave electron-phonon (EP) pairing. These exam- 
ples suggest there are important aspects of two di- 
mensionality (2D) for conventional superconductors 
that are yet to be understood. 

A close analogy for MgB2, structurally, electron- 
ically, and regarding superconductivity, is graphite. 
Graphite has the same C layer structure as B has 
in MgB2. The layer stacking that occurs in graphite 
is central to its semimetallic character but is not 
important in our analogy. Graphite is isoelec- 
tronic with MgB2 ; previous studies have established 
that the Mg atom is effectively ionized. Finally, 
graphite becomes superconducting up to 5 K, but 
only when doped (intercalated). || Both graphite 
and MgB2 have planar sp 2 bonding, and in graphite 



it is well established that three of carbon's four va- 
lence electrons are tied up in strong a bonds lying in 
the graphite plane, and the other electron lies in non- 
bonding (p z ) 7T states. Electron doping of graphite, 
achieved by intercalating alkali atoms between the 
layers (Na is the most straightforward case) leads to 
occupation of otherwise unfilled ir states and leads 
to T c as high as 5 K. MgB2, on the other hand, 
superconducts at ^40 K at stoichiometry. 

The light masses in MgB2 enhance the phonon 
frequency (ujph oc Af^ 1 / 2 ) that sets the tempera- 
ture scale of T c in BCS theory. Even considering 
this tendency, there must be some specific feature(s) 
that produces such a remarkable T c , and moreover 
does so with no d electrons, nor even the benefit of 
a density of states (DOS) peak. [|l0| , |ll"f In this paper 
we identify these features: 

(1) hole doping of the covalent a bands, achieved 
through the ionic, layered character of MgB2, 

(2) 2D character of the a band density of states, 
making small doping concentrations give large 
effects N(e^)^ m* /nh 2 independent of rih, and 

(3) an ultrastrong deformation potential of the a 
bands from the bond stretching modes. 

In this paper we exploit the similarities between 
graphite and MgB2 to build a basic understanding 
of its electronic structure, then focus on the differ- 
ences that are connected to strong EP coupling, and 
finally provide an estimate of T c that indicates the 
picture we build can account for the observations. 
This picture further predicts that electron doping 
by ~0.2 carriers per cell will strongly affect T c ad- 
versely, a result that has been reported by Slusky 
et al. |1J] Although Hirsch has also focussed on the 
hole character of the a bands, [jl3| his emphasis is 
otherwise quite different from that described here. 

Calculations of the electronic structure have been 
done using the linearized augmented plane wave 
(LAPW) method Q that utilizes a fully general 
shape of density and potential, as implemented in 
the WIEN97 code. [fi~5| Experimental lattice con- 
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stants of a=3.083 A, c=3.521 A were used. LAPW 
sphere radii (R) of 2.00 a.u. and 1.65 a.u. were cho- 
sen for the Mg and B atoms, respectively, with cut- 
off RK maa; =8.0, providing basis sets with more than 
1350 functions per primitive cell. The generalized 
gradient approximation exchange-correlation func- 
tional of Perdew et al. Jl6| was used in the present 
work. 

(a) 



4.0 













E F 












r m k r a l 












E F 














r m k r a l 



FIG. 1. Band structure along main hexagonal sym- 
metry lines, for (top) MgB 2 , (middle) □ 2+ B 2 , (bottom) 
primitive graphite C2. The planar a states, highlighted 
with larger symbols, fall in energy in this progression, 
and only in MgB 2 are they partially unoccupied. The 
point A = (0,0,7r/c) is perpendicular to the (k x ,k y ) 
plane. 

The band structure of MgB2 is shown in Fig. |l| 
(top panel) in comparison with that of primitive 
graphite (bottom panel) with a single layer per cell 
like the B2 sublattice in MgB2. For each two distinct 
sets of bands are identifiable: the highlighted sp 2 (a) 
states, and the p z (tt) states. The striking difference 
is in the position of the a bands, which is evident 
in Fig. [l| Whereas the a bonding states are com- 
pletely filled in graphite and provide the strong co- 



valcnt bonding, in MgB2 they are unfilled and hence 
metallic, with a concentration of 0.067 holes/B atom 
in two fluted cylinders surrounding the T-A line of 
the Brillouin zone. |l0| There are correspondingly 
more electron carriers in the tt bands. This decrease 
in occupation on the strongly bonding a bands par- 
tially accounts for the greatly increased planar lat- 
tice constant of MgB2 (3.08 A) compared to graphite 
(2.46 A). Our results agree with previous conclusions 
that MgB2 can be well characterized by the ionic 
form Mg 2+ (B 2 ) 2 -. 

To identify the origin of the relative shift of the 
a and tt bands by ~3.5 eV between graphite and 
MgB2, we have considered a fictitious system D 2+ B2 
in which the Mg ion is removed but the two electrons 
it contributes arc left behind (and compensated by 
a uniform background charge). The band structure, 
shown in the middle panel of Fig. [l], is very sim- 
ilar, except the energy shift of ~1.5 eV downward 
with respect to MgB2 completely fills the a bands, 
as in graphite. This shift is the result of the lack 
of the attractive Mg 2+ potential in MgB2, which is 
felt more strongly by the tt electrons than by the in- 
plane a electrons: the attractive potential of Mg 2+ 
between B2 layers lowers the tt bands, resulting in 
a — > tt charge transfer that drives the hole doping of 
the a bands. Belashchenko et al. Jjj]] have also con- 
sidered a sequence of materials to come to related 
conclusions about the band structure, but they did 
not use isoelectronic systems as has been done here. 

The a bands are strongly 2D (there is very little 
dispersion along T-A), but it will be important to 
establish the magnitude and effects of interplanar 
coupling. The light hole and heavy hole a bands in 
MgB2 can be modeled realistically in the region of 
interest (near and above £p) with dispersion of the 
form 

kg. ~\~ k^. 

£ k = £o - - 2t±cos(k z c), (1) 

Zm* 

where the planar effective mass m* is taken to be 
positive and t± = 92 meV is the small dispersion 
perpendicular to the layers. The light and heavy 
hole masses are m^ h /m=0.20, m* hh /m= 0.53, and 
the mean band edge is e o =0.6 eV. In general, the in- 
plane (v xy ) and perpendicular (v z ) Fermi velocities 
are expected to be anisotropic: v xy ~ kp/m*, v z ~ 
2ct± where t± is small. Near the band edge (kp < 
2m*ct±) this anisotropy becomes small, and this is 
roughly the case in MgB_2. The tt bands are also 
effectively isotropic. [ |10|Jllf| 

Now we discuss why the quasi- 2D character of 
the a bands is an important feature of MgB2 and 
its superconductivity. Neglecting the k z dispersion, 
the 2D hole density of states is constant: N£(e) 

— m ' h J ^i hh — 0.25 states/eV-cell, independently of 
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the fact that the hole doping level is small. The k z 
dispersion has only the small effect displayed in Fig. 
||, where the discontinuity in the quasi-2D DOS is 
seen to be broadened by ~ 2t±. For MgB2 the a 
band contribution to N(ejr) is reduced by about 10% 
by k z dispersion. 
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FIG. 2. Density of states (solid lines) of the light hole 
band, the heavy hole band, and the total, for the model 
of Eq. (1). Note that N(e) drops rapidly only 0.4 eV 
above ef- The dashed lines give the 2D (t±=0) analogs. 

If superconductivity is primarily due to the ex- 
istence of holes in the a band, and we provide evi- 
dence for such a picture below, then the DOS in Fig. 
H suggests that electron doping will decrease N(ef). 
The decrease will be smooth to a doping level corre- 
sponding to an increase by 0.4 eV of the Fermi level, 
i Then N(e.f) should drop precipitously with further 
doping. A rigid band estimate gives a value of n^ ;Cr 
=0.08 electrons necessary to fill the a bands to the 
crossover region. Electrons must also be added to 
the 7T bands that lie in the same temperature range, 
making the total doping level n cr w 0.25 per cell. 

Fig. H therefore suggests a critical region of elec- 
tron doping around 0.25 carriers per cell. Slusky et 
al. have reported a study of such electron doping 
in the Mgi_ a; Al a ;B2 system. T c drops smoothly up 
to .x=0.1, beyond which point a two phase mixture 
of B-rich and Al-rich phases occurs. At x — 0.25 
and beyond, a single phase, non-superconducting 
system is restored, strongly supporting our picture 
that filling the a bands will destroy superconductiv- 
ity. Using the rigid band picture, filling them de- 
creases T c moderately initially, but as the a bands 
become nearly filled, the coupling decreases abruptly 
and T c vanishes, as observed. Although our results 
do not bear directly on the two-phase question, we 
note that this occurs just as the a bands are filling. 
At this point N(ef) is dropping, which normally fa- 
vors stability. The observed instability suggests that 
a very small density of a holes that are very strongly 
coupled to the lattice is what underlies the lattice in- 
stability. 

Now we address the question of EP coupling 



strength A. The calculated value of N(e^) =0.71 
states/eV-cell correponds to a bare specific heat co- 
efficient 70 = 1.7 mJ/mole-K 2 . The available ex- 
perimental estimate is Q Jexp = (1 + ^)"fo — 3 ± 
1 mJ/mole-K 2 , giving X exp ~ 0.75. Kortus et al. 
have used the rigid muffin tin approximation |lS| ] 
to obtain an idea of the coupling strength. How- 
ever, this approximation is not well justified in sp 
metals and neglects distinctions between bands of 
different character that we expect to be crucial for 
high T c . If one assumes the wavevector dependence 
of coupling is not strong, there is another simple way 
to identify strong coupling using deformation poten- 
tials V = Ae^/AQ due to frozen-in phonon modes 
with mode amplitude Q. |19| The underlying con- 
cept is that a phonon that is strongly coupled to 
Fermi surface states will produce a large shift in Ek 
for states near the Fermi level. |l9) 

We have studied these deformation potentials for 
the A;=0 phonons Bi g , E2 ff , A.2u> Ei u whose energies, 
86, 58, 48, 40 meV respectively, have been calculated 
by Kortus et al. |l(]] This mode (and others) destroy 
the symmetry of the crystal. The Brillouin zone re- 
mains unchanged, however, and we plot the bands 
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along the same directions as in Fig. hi This mode, 
as well as the Bi g mode, involves only out-of-phase 
motion of the two B atoms in the primitive cell, E2 9 
involving in-plane vibrations (bond stretching), Bi s 
having displacements perpendicular to the plane. In 
Fig. the actual rms B displacement of Aub = 
0.057 A was used, to provide a clear picture of the 
very large effect of EP coupling strength. 




FIG. 3. Band structure with frozen-in F<2 g mode, plot- 
ted along the same lines as in Fig. 1 to facilitate com- 
parison. (Note that the point labelled K is no longer a 
symmetry point.) The E2 3 phonon breaks the symmetry 
and splits the a bands by 1.5 eV along T-A. 

The Fi2 g phonon strongly splits the a band nearly 
uniformly all along the T-A line and near ef, with 
the "gap" opening A£ gap /Au B = 26 eV/A= 2V E2g ; 
note that it will be the square of T> that enters the 
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expression (below) for A. In stark contrast, the Bi g , 
A2 U , and Ei„ modes produce no visible effect in the 
bands; we estimate that their deformation potentials 
are at least a factor of 25 smaller. Thus the only 
significant deformation potential is for the a bands 
(-7T band shifts are always small), and only due to 
the E2 S mode. This one is extremely large, suggest- 
ing that non-linear coupling may even be occurring. 
This strong coupling should be observable as a large 
linewidth, and the superconducting gap can be 
expected to be larger on the a Fermi surface sheets 
than on the tt sheets. 

To estimate the coupling from this mode alone, 
we use Eq. (2.34) of Kahn and Allen jl9| for the 
EP matrix element M in terms of T>E 2g = 13 eV/A, 
including only the two a bands N^ = 0.11 eV _1 per 
spin (from Fig. 2) and the two E2 S modes: 



^E 2g =N h (e F ) 



v=l,2 



\MU 2 



(2) 



= 2N h (s F )[ 



2M B u 2 



1.0. (3) 



.3=1,2 



The sum on j runs over the two moving B atoms 
of the E2 g mode. Using this value of A, the E2 9 
mode frequency, and the Allen-Dynes [^0| or McMil- 
lan equation with Coulomb pseudopotential fi* = 
0.10 - 0.15, the resulting range is T c = 32-46 K. 
There is considerable uncertainly in this estimate, 
but it seems quite plausible that coupling of the a 
band to the bond stretching mode may provide most 
of the coupling to account for the observed value of 
T c . 

Now we summarize, beginning with the effects of 
low dimensionality mentioned in the introduction: 
(1) the B2 layers provide a strong differentiation be- 
tween B states (a vs. it) that results in the Mg 2+ 
layer giving a 3.5 eV a — it energy shift, driving self- 
doping of the a bands, (2) due to their 2D dispersion, 
the contributions of the a bands to N(ef) is almost 
independent of the doping level; specifically, for low 
doping level N— h(ep) is not small. Then, the strong 
covalent nature of the a bands leads to an extremely 
large deformation potential for the bond stretching 
modes. Judging from the very small deformation po- 
tentials of the other three T point modes, the bond- 
stretching modes dominate the coupling. 

W.E.P. is grateful to P. C. Canficld for conversa- 
tions and for early communication of manuscripts. 
This work was supported by Office of Naval Research 
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Note added: Since this paper was submitted, two 
reports of EP coupling in MgB2 that support our 
results [21 22jl have come to our attention. 
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